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CONTROLLABILITY OF NETWORKS OF ONE-DIMENSIONAL SECOND 
ORDER P.D.E. AN ALGEBRAIC APPROACH 

FRANK WOITTENNEK AND HUGUES MOUNIER 

Abstract. We discuss controllability of systems that are initially given by boundary coupled 
p.d.e. of second order. Those systems may be described by modules over a certain subring 1Z 
of the ring A4o of Mikusihski operators with compact support. We show that the ring TZ is 

QQ ' a Bezout domain. This property is utilized in order to derive algebraic and trajectory related 

f - *j ' controllability results. 
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1. Introduction 

The solution of control design problems is, in general, preceded by a controllability analysis of 
the system under consideration. While for linear finite dimensional systems, both algebraic and 
analytic controllability notions are used in parallel, the analysis of infinite dimensional systems is 
dominated by (functional) analytic methods [9]. The latter approach has proven to be useful, in 
particular for the analysis of state space controllability, i.e., the possibility of steering the system 
under consideration from a given initial state to a desired final state. For example, controllability 
of the same class of systems as considered in the present contribution has been analyzed this way 
_C ■ in [10,20]. However, when focusing on trajectory tracking problems, the behavioural controllability 

notion due to Willems [45] is an interesting alternative to classical state space controllability. The 
connections between this approach and the algebraic system properties have been pointed out in the 
past for several classes of distributed parameter systems, as for so called multidimensional systems 
[33,49], for delay systems [17,30], and more general convolutional systems [43,44]. In addition, 
the algebraic viewpoint is closely related to parametrization of trajectories of the system under 
consideration. This constructive nature makes the approach very attractive for applications, in 
particular for open loop control design. Finally, paying attention to particular structural properties 
of the models under consideration may result in a deeper understanding of the respective class of 
systems. 

From the algebraic (module theoretic) viewpoint, as used within this contribution, a linear 
system is a finitely presented module. This notion was first introduced by Fliess for linear finite 
dimensional systems in [11]. For this class of systems, the freeness of the module corresponds to the 
flatness of the system under consideration in the sense of the theory of nonlinear finite dimensional 
systems while its basis corresponds to a flat output [12]. Moreover, torsion freeness, i.e., the absence 
of autonomous subsystems, is equivalent to freeness. The module theoretic approach is applicable 
to systems with distributed parameters and lumped controls as well: The convolutional equations 
associated with a given boundary value problem serve as defining relations for the system module, 
the latter defined over a certain ring of generalized functions. Such coefficient rings are generally not 
principal ideal domains. For this reason, the two basic controllability related properties, torsion 
freeness and freeness, are not necessarily equivalent. An approach to circumvent the problems 
caused by this "lack of structure" is the concept of 7r-freeness, which relies on localization and was 
at first developed for linear delay systems [13]. This way a basis can be introduced at least within 
an appropriate extension of the module under consideration. The approach has been proven to 
be very useful for both trajectory planning and open loop control design [21,31,36,37,39,46,47]. 
Nevertheless it seems to be difficult to compare such purely algebraic controllability notions to the 
behavioural ones. For this reason, within the present contribution we do not use localization. 

This paper addresses the development of an algebraic approach to the controllability of networks 
of spatially one-dimensional parabolic and hyperbolic constant coefficient p.d.e. of second order. 
Here, by a network we understand a system consisting of several branches each of which is governed 
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2 F. WOITTENNEK AND H. MOUNIER 

by a system of p.d.e. and which are coupled via the boundary conditions. We investigate, through 
algebraic properties of the coefficient rings obtained for considered class of p.d.e., the related 
controllabilities of the associated system module and establish some controllability results including 
module theoretic and behavioral ones. In accordance with [38,39,46], we use the general solution 
of the Cauchy problem with respect to space in order to rewrite the given model as a linear 
system of convolutional equations. The latter are regarded as the defining relations of a finitely 
presented module. It turns out that the coefficient ring of this module, a subring of the ring of 
Mikusiiiski operators with compact support introduced in [1], is a Bezout domain, i.e., every finitely 
generated ideal is principal. An algorithm enabling us to calculate the generator of a given finitely 
generated ideal is presented within this paper. This latter result is strongly inspired by those 
derived in [2,17] for particular rings of distributed delay operators which in our setting may arise 
from the wave equation. The derived properties of the coefficient ring allow us to decompose the 
system module into a free module and a torsion module. Finally, from these algebraic results, we 
deduce the trajectorian controllability of the free submodule in the sense of [15] and its behavioural 
controllability in the sense of [45]. 

The paper is organized as follows. In Section [2] we introduce the class of models considered as 
systems of p.d.e. which are coupled via their boundary conditions. We show how to pass from this 
model to a system of convolution equations giving rise to our module theoretic setting. Section 
|3] is devoted to the study of the coefficient ring of this module. In Section [4] we obtain several 
controllability results for the systems under consideration. Finally, in Section [5l we apply the 
method to a system example of two boundary coupled p.d.e. 

2. Boundary value problems as convolutional systems 

2.1. Models considered. We assume, that the model equations for the distributed variables in 
iPi, . . . , wi and the lumped variables in u = (u\, . . . , u m ) are given by 

d x w l — A l w i + BiU, wi : fi$ — > B 2 , u G B m 

At G (M[s]) 2x2 , B t G (M[s]) 2xm 

where B denotes an appropriate space of Boehmian4j {see, e.g., [6, 25, 27] and App. |SJ and s 
is the differentiation operator with respect to time. The assumptions which are crucial for the 
applicability of our approach are twofold. First, we assume that all the matrices A\, . . . , Ai give 
raise to the same characteristic polynomial, namely, 

(lb) det(AI- A,) = A 2 -a, a = as 2 + bs + c ^ 0, a, b, c G M, a ^ 0. 

Additionally, we require the intervals Qi, . . . , $7/ of definition of the above differential equations to 
be rationally dependent. More precisely, we assume the fij (i = 1, . . . , I) to be given by an open 
neighbourhood of 

(lc) Qi = [xi,o,Xi t i], h = Xj,i - Xifi = qd, % G Q, i G R. 

In the following, and without further loss of generality, we assume Xi t o = 0. The model is completed 
by the boundary conditions 

(Id) Yl L i™M + RiWi(£i) + Du = 



7=0 



<jx2 



where D G (M[s])« xm and U,R % G (R[s)) 

Remark 2.1. In a more general setting, instead of the boundary conditions (|ld|). one could con- 
sider auxiliary conditions of the form 



I 

E 

i=0 



Q t (wi) + Du = Q. 



Alternatively, one could use other spaces of generalized functions given in the inclusion chain £ p C £ C T>' C 
T>' C B. Here, as usual, £ denotes the space of infinitely differentiable functions, T>' are the Schwartz distributions. 
Moreover, £ p and T>' are spaces of Gevrey- functions and Gevrey ultra-distributions of order p, respectively. However, 
in order to avoid distinctions of cases, we shall use the space of Boehmians. 
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Here, 

v H „ 

Qi{Wi) = ^ L i,3 W i( a ijt) + ^2 / Qi,j(x)lVi(x)dx 
3=0 3=1 Jn '-i 

withL itJ e (K[s])? x2 , Q* 3 e (R[s,a;]) gx2 , ft 4 d fi^ - [A,j,i^, ft,j,2^], Oi^faj, e Qnft,, ftV eN. 

2.2. Solution of the inital value problem. This section deals with the solution of a single 
initial value problem of the form (|la[) with initial conditions given at x = £, i.e., 

(2) d x w = Aw + Bu, w(£) = w^ 

with A, B having the same properties as A4, Bi (i = l,...,l) introduced within the previous 
section. To this end, we start with the initial value problem 

{d 2 x - a)S(x) = 0, S(x) = 0, (d x S)(x) = 1, 

associated with the characteristic equation (jlbj) . It is well known that this equation has a unique 
operational solution as long as the principal part d 2 — as 2 of the differential operator d 2 — a is 
hyperbolic w.r.t. the parallels of x — 0. This was implicitly required above by assuming a Js 
in (|lb[) . Moreover, under these assumptions the operator S as well as its derivative C = d x S 
correspond to infinitely differentiable functions mapping O to the ring M of Mikusihski operators 
with compact supporio (see App.lBland [1] for results related to the support of Mikusihski operators 
and App.lAlfor explicit expressions for C{x) and S(x)). 

Using the above defined operational functions one easily verifies that the (unique) solution 
x <— > $(x, £) of the initial value problem 

d x $(x,0 = A$(x,£), $(0 = /, 
with / denoting the identity, is given by 

(3) $(x,S) = AS(x-t)+IC(x-S). 
From the uniqueness of the solution one deduces the addition formula 

For A the companion matrix of the characteristic polynomial, i.e., 

U) A -(° A $( X £)-( G W S W 

(4) A - [a 0)> *&>&- [aS(x) C(x] 

this yields in particular 

(5) C(x + y) = C(x)C(y)+aS(x)S(y), S(x + y) = C(x)S(y) + S(x)C(y). 
The solution of the initial value problem associated with the inhomogeneous equation 

(6) d x *{x,S)=A*(x,Z) + B 

with homogeneous initial conditions, prescribed at x = £, is obtained using the well known variation 
of constants method. This yields 



(7) *(z,0=/ $(x,()d{B. 

Thus, the general solution of the initial value problem @ reads 

w(x) = $(x, £)W£ + $f(x, £)u. 

The entries of the matrix $ belong clearly to C[s, C, S]. Contrary, according to 0, the entries 
of \P may contain also the integrals of S and C. However, if o ^ then A is invertible over C(s). 
Thus, using the fact that d z Q(x, z) — — $>(x, z)A those integrals can be expressed as 

f $(x, C)dC = - f d c <P(x, QA-'dC = ($(», - I) A- 1 . 



More specifically, instead of stating that S(x), C(x) £ .Mo, one could distinguish the cases a > and a = 0. In 
the first case both operators, S(x) and C(x), correspond to distributions of order zero with compact support, while 
in the latter case they correspond to ultra-distributions of Gevrey order 2 and support in 0. Both of these spaces 
may be embedded in .Mo- 
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Choosing A as in ^ one obtains in particular 

C{Qdx = S(x), I S{()dx = (C(x) - l)/cr. 
Jo 

Later on, the latter equations will essentially ease our controllability analysis. 

2.3. A module presented by a system of convolution equations. In the previous section 
we have discussed the solutions of the initial value problems associated with the equations (jlaj) . In 
the sequel these results are used in order to define an algebraic structure representing the model 
under consideration. 

Substituting the general solutions of the initial value problems into the boundary conditions, 
one obtains the following linear system of equations: 

(8) w{x) = W^{x)c t , P € c € = 0. 

Here, | = (fc, . . -,£»), & = «(fc) . . . w f^), „T )) 



'*i(s,&) *i(x,& 

'•• : 



w*=l ••. o i h n = (^,i---^+i) 



with 



p w+1 = r> + ]r l 4 * 4 (o, eo + R&iiii, &). 

A suitable algebraic object for the representation of the model under consideration should 
contain all the system variables, i.e., the distributed variables in w, their values as well as their 
(spatial and temporal) derivatives. Moreover, it should reflect not only the structure imposed by 
the original boundary value problem (|la[) but also that imposed by the solution of the according 
initial value problems, i.e., by equation l[5|). In order to analyze the model we will, therefore, use a 
module generated by the variables collected in c^ with the presentation given in (JSj> [13,14,16,29]. 
The choice of the coefficient ring, which has to contain at least the entries of W$ and Pg, is discussed 
below. 

According to the previous section, the entries of W$ and those of P$ are composed of the 
functions C and S mapping ItoA^o and all the values of these functions. Moreover, the matrices 
may involve values of the spatial integrals of C and S, too. Thus a possible choice for the coefficient 
ring is the ring K^[s, 6, 6 1 ]. Here, for any X C R, ft£ = [6 X , &{}, with 

e = {C,S}, 6 x = {C(z£),S(z£)\zeX}, 

e 1 = {c 1 , s 1 }, 6| = {c'izi), sV)l* e x}, 



S J (x) = / 5(C)dC, C 7 (x) = / C(CR 



defined as in (fTc|) . and 



Inspired by the results given in [2, 17,30], and in view of the simplification of the analysis of 
the module properties, instead of the ring 1Z^, we will use a slightly larger ring, given by 72.R = 
C(s)[6 R }nM . 

Definition 2.1. The convolutional system E associated with the boundary value problem ([T]) is 
the module generated by the elements of c^ over 1Z = 7£r[<5, & 1 ] with presentation matrix P^. By 
Er (resp. Eq,) we denote the same system but viewed as a module over 72-r (resp. TZq). 

One easily verifies that E does not depend on the choice of £ (cf. [39,46]). In view of the assumed 
mutual rational dependence of the lengths £\, . . . ,£i, for the analysis of the system properties, we 
will start with the system Eq, i.e., a system containing the values of the distributed variables at 
rational multiples of £ only. However, having analyzed the properties of Sq, we may pass to Er 
(resp. E) by an extension of scalars, i.e., Er = 7£r ®Kq Eq (resp. E = 1Z (3iz Q Eq). 
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3. The operator ring IZq is a Bezout domain 

In this section we study the structures of the ideals within the rings IZq. To this end, we first 
establish some results on the ideals in C(s)[©q] and C(s)[©n]- 

3.1. Ideals in C(s)[6q] and C(s)[Sn]- This section is devoted to the analysis of the ideals in 
C(,s)[6q] and C(s)[Sn]- The obtained results will be the key for the analysis of the properties of 
TZq and 72.pj introduced in Section 12.31 

In the following, we will replace C(s) by any field k. Moreover, the generating set 6x may be 
replaced by any set G5x the elements of which satisfy the addition formulas derived in Section 12.21 

Definition 3.1. Let k be a field and X an additive subgroup of R. By ®x we denote the set 
{S a , C a \a € X} the elements of which are subject to the following relations (a G k): 

(9a) C a Cb ± aS a Sb = C a ±b, S a Cb ± C a Sb — S a ±b 

and 

(9b) Co = 1, S = 0. 

From the above definition one easily deduces 

(10a) G a = C_ a , o a = — 0- a 

as well as 

(10b) 2C a Cb — C a +b + C a -b, 2aS a Sb — C a +b — C a -b, 2C a Sb = S a +b — S a -b- 

Moreover, any element r £ fc[(5x] can be written in the form 

n 

(11) r = y^ a, ai C ai + b ai S ai , neN, a ai ,b ai € k, a;€X + 

i=0 

where X + = {\a\ : a <G X}. Finally the units in fc[©x] belong to k. 

In the following, it is necessary to distinguish the cases where the equation A 2 — a = has a 
solution over k or not. For our application this is clearly equivalent to the question whether the 
roots of the characteristic equation (|lb|) belong to K[s]. The necessity to distinguish these cases 
is explained by the following simple example which, in addition, shows that the cases X = N and 
X = Q need to be analyzed separately. 

Example 3.1. Consider the ideal J = (a, b), a — Si, b = C\ + 1. Over fc[6tj] we have 

a = Si = 2Ci/ 2 Si/ 2 , b = C\ + 1 = 2Cjy 2 . 

Thus, both generators belong to (Ci/ 2 ) which, conversely, belongs to 3 since 2Ci/ 2 — —o-S\i 2 a + 
Ci/ 2 b. The ideal 3 is, therefore, generated by C\/ 2 which does not belong to fc[6jj] */ -^ 2 — °~ * s 
irreducible over k. However, the situation is different if ■s/a belongs to k. From the relations given 
in (0) and U0\) . it follows immediately that 

(Ci/2 + y/a S1/2 )(Ci/ 2 - y/a Si/ 2 ) = 1. 
Over k[&o\, C\/ 2 can be factorized as 

C 1/2 = (C 1/2 + V^S 1/2 )(l + a - ^Si)/2. 
The element C\/ 2 is, thus, associated with 1 + C\ — \f5S\ which indeed belongs to /c[6n]- 
3.1.1. The polynomial A 2 — a is reducible over k. 
Proposition 3.1. The ring &[<5n] is a PID. 

Proof. From the addition formulas given in ([9]) and (fT0|) . it follows that A:[©n] is isomorphic to 
k[S±, C{\ which, in turn, is isomorphic to k[z _1 , z] by 

7- 1 - z 1 
61 1-^ r , Gi^Z + Z . 

A 
The latter ring is Euclidean with the norm function given by the difference of the degrees of the 
monomials of maximal and minimal degree w.r.t. z. □ 

Corollary 3.1. The ring fc[(8<j] is a Bezout domain. 
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Remark 3.1. Note that, for the same reason as given in Remark \3.3l fc[<S(Q>] is not a PID. 



Remark 3.2. Having in mind that in our application a is given according to equation Ub\) and 
k = C(s), for y/a G k, say yfa = A, the operators C(x) and S(x) introduced in Section \2.2\ arc 
constructed from point delays: From A = Sy/a — a, a G C we obtain 2C(x) — e sv ^~ Q + c-( s Va- a ) , 
Note that in this case our results are simply a restatement of those presented in [2, 17]. 

3.1.2. The polynomial A 2 — a is irreducible over k. As indicated by Example 13. 1[ the second case, 
i.e., the equation A 2 — a = has no solutions over k, is much more challenging than the first 
one. There, the ring A;[<5n] corresponds basically to the ring Q[x, y]/[x 2 + y 2 — 1] of trigonometric 
polynomials which is obtained from fc[©pi] for a — — 1 and k = Q. The latter ring is lacking the 
pleasing properties of a PID or even a Bezout domairo However, Example 13. II suggests, that the 
difficulties can be circumvented when allowing to halve the argument, i.e., working with &[(8q] 
instead of fc[©N]- 

Definition 3.2. For any nonzero r G fc[<3x] the norm v(r) is defined as the highest a G X + such 
that at least one of the coefficients a a and b a in (jllj) is nonzero. 

Lemma 3.1. Let S the multiplicative subset ofk[&fi] consisting of all the elements such that either 
the coefficients with odd or those with even indices vanish. More precisely, any element s of S can 
be written as 

s = J2 a Sti Ci + b a ,iSi, I s = \v(s) -2»|ieZ, 0<*<^pi 

Letp,q G S the norms of which are strictly positive. Without loss of generality assume v(p) ^ K?)- 
Consider the ideal J — (p, q) generated by p and q. Then there exists p,q G S with 3 — {p, q) and 
either v(p) > v{p) ^ v(q) or q = 0. 

Proof. In the following, three different cases are considered. 

Case 1. If v(p) > v(q) one can apply a division step similar to that of polynomials. More 
precisely, we will show that there exists r,h G S with either r = or v{r) < i/(p) such that 
p = qh + r. Then we may set p = q, q — r (or vice versa) to complete the discussion of the first 
case. 

In order to show that r, h with the claimed properties exist set 

h = a h CA + bh.SA, A = v{p) - v(q) 
where the coefficients ah , bh G k have to be determined appropriately. It follows 



s = hq = 22 (( a ha q ,iCiCA + bhaq^C'iSA) + (a h b q ^SiC'A + &/i^,jS'j<S'a) j 
— 7T- J^ ( {o-a h a q ,i + b h b q ,i)CA+i + (cra h a q ^ - b h b q ^)C a-i 



iei„ 



+ o-(b h a q .i + dhb q ,i)SA+i + o~{bha q ,i - a h b q ^)SA-i\ 






a s,v{p) = 7rV ja hO'q,v(q) +W^(j))i ^s,u{p) = ^ \phO>q,v{q) + °>hfcq,v{q)) ■ 



/ Ojs,iCi + St ibi 

where the leading coefficients are given by 

— (o-aha q , v ( g ) + bhb qil/ (q)) , & s ,i/(p) = ~l 

From this equation and from r = hq — p the norm of r is smaller than that of p if and only if ah, 
bh satisfy 

( 12 ) ( a q,v(q) a ~%,»(<iA f a A _ 2 ( a P^(p) 

\bqMq) a q ,u(q) J \b h J \ b p^(p) 

By the definition of the norm at least one of the coefficients a q , v (q) an d a q,v(q) is nonzero. Since 

i 2 - h 2 



additionally, yfa G" k it follows aa 2 v , % — b 2 v , •. ^ and an, bh can be alway chosen according to 



Actually, the trigonometric ring is a Dcdckind domain [7,8,32]. 
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Case 2. If v(p) = u(q) and for some c 6 k the equations a q , u ( q ) — ca pu ^, b qv ( q ) — cb PjV ^ 
hold, the ideal 3 is generated by p — p, q = q — cp where v(q) < v{p)- If q = the proof is complete 
otherwise we can proceed according to the first case with the pair p, q instead of p, q. 

Case 3. If v(p) = v{q) but we are not in the second case set 

(13a) (p q) T =A 1 (p q) T 

(13b) (p q) T =A 2 (p qf 



__ I %),n Op,n \ ,, I (-'l 0~t>\ 



with 

Obviously, p, q belong to the ideal generated by p, q. Both matrices, Ai and A%, are invertible, the 
first one since otherwise we would be in the second case, the latter one since, by © , its determinant 
equals 1. Thus, (p,q) = (p,q) = (p,q). 

It remains to show that the norms of p and q are both smaller than n. From equation Q13ap one 
obtains v{p) = v{q) = n with ap tV ( q ) = b q . u ( q ) — 1, &p,i/(g) = a q,v(q) — 0- From (|13b|) one has 

p = C\C n - aSiS n + ^2 Ci (ap^Ci + bpjSi) - aSi (a q ,id + b q ,iSi) 



q — C\S n — S\C n + 2_^ Ci(a q iCi + b q .iSi) — S\{a,p^Ci + bp-jSi) 



E 



with I* = I p \{n}. The norms of the sums in the above expression are at most n — 1 while for the 
leading terms one obtains according to (p?aj) 



WW — <Jd\b n — C n ~\, C\S n — S\C n — bn—l- 

Thus, the norms of p, q cannot exceed n — 1. □ 

Lemma 3.2. Letp,q eSc fc[®NJ with v{p) Js v(q). Then there exists p,q e Sfl (p, q) such that 
(P, q) = (P, «) a«rf K<?) < K«); ^(p) < K«) or g = 0. 

Proo/. By Lemma 1531 (p. q) = (p*,q*) with v(p) > v{p*) ^ v{q*) or g* = 0. In the latter case the 
claim has been proved. Otherwise, repeat the above argument p*,q* until we are in the claimed 
situation which happens after at most v{p) — v(q) + 1 steps. □ 

Proposition 3.2. Any ideal 3 in k[<3ff] generated by a subset of S is principal. 

Proof. Step 1. We show that up to multiplication with units there is only one element q of lowest 
norm v{q) = n in 5* fl 3. To this end, assume there are at least two such elements, say p and q. By 
Lemma 13.21 there exist p,q € S with (ft, q) — (p, q) where n > v{p) ^ v(q) or n ^ v{p) and q = 0. 
Since n is the lowest possible norm for an element of 3 ("I S, only the case n = v(ft) and q = 
remains. But this can happen only if we are in case 2 of Lemma 13.11 having p = p and q = cp, 
c£ k y . 

Step 2. We now show that any element of (5 belongs to (q) where q is defined as in the first 
step. To this end chose any element p from ©. Applying case 1 of Lemma 13.11 several times one 
gets p — hq + r, v(r) ^ n, r e S. Since, by assumption, q has the smallest possible norm, it 
follows v(r) — n or r = 0. This in turn yields r — cq, c 6 k according to Step 1. Finally, we have 
p = (h + c)q and, therefore, 3 = (q). □ 

Proposition 3.3. Any finitely generated ideal 3 in A;[(5q] is principal, i.e., fc[0Q] is a Bezout 
domain. 

Proof. Let 3 = (ri, . . . , r m ) for some m G N. Write the generators according to (jTTJ), i.e., 
(14) Tj =^2a aiti C a:jii +b ajti S ajii , Uj € N, a aitj ,b aiJ € k, a iyj eQ + . 

i=0 

Let <i be a common denominator of all the ckjj occurring in these equations. Then the generators 
of 3 can be identified with elements of the subset S defined in Lemma (|3.1[) of the Ring R® via 
the embedding E : i?pj — > i?Q which is defined by C2 1— ► Ci/^i ^ | -^ S'l/d- Let f\, . . . , r m elements 
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of Rn the images of which are r\ , . . . , r m . The ideal 3 generated by f i , . . . , f m is principal by 
Proposition 13. 21 Consequently, 3 is generated by the image of the generator of 3 under E. □ 

Remark 3.3. Note that neither A;[®q] nor /c[©n] are principal ideal domains (PID). The first is 
not Noetherian: As an example for an ideal that is not finitely generated take {{Si/ 2 n\n £ N}). 
Moreover, fc[65jj] * s n °t a PIT) since there are finitely generated ideals that cannot be generated by 
one single element: The ideal (Si,Ci + 1) viewed as an element of fe[(8{j] is generated by C\i 2 
which does not belong to fc[(5pj]- 

3.2. TZq is a Bezout domain. We are now in position to prove that TZq is a Bezout domain. 
After the preparation done in the previous subsection the remaining steps are very similar to those 
given in [2, 17]. In particular, the proof of Lemma 13.31 which prepares Theorem 13.11 is strongly 
inspired by [2, Theorem 1]. 

Lemma 3.3. For two coprime elements p,q € TZq there exist a,b € TZq such that ap + bq = 1. 

Proof. By Prop. l3~3"l (resp. Cor. l3.ip C(s)[6q] is a Bezout domain. Thus, there exist a, b £ C[s, 6] C 
TZq such that ap + bq = h where h € C[s]. Write h as product h — Y[i=i( s ~ s i) an d proceed by 
induction (we do not assume s, ^ Sj). 

Assume there exist a, 6 € TZq with ap + bq = Ili=i( s — s «)- I n the following, for any 7 € A^o we 
set 7 = = 2 , (7)(sat), with the entire function ^(7) being the Laplace transform of 7 given according 
to App. IB. 21 By the coprimeness of p, q and by Lemma IB. 31 in the Appendix, a* and 6*, defined 
by 

(s-s N )\q \— -, {s-S N )\p 

, b*=\ P( s ~ s ^) 
(s~s N )\q . f s - SN )\ p 

v S — Sn 

belong to TZq. One easily verifies that pa* + qb* = Y\ i=1 (s — Si). Applying this step N times 
completes the proof. □ 

Theorem 3.1. The ring TZq is a Bezout domain, i.e., any finitely generated ideal is principal. 

Proof. We show that any two elements p, q € TZq possess a common divisor c that can be written 
as linear combination of p, q. (It is then the unique greatest common divisor of p and q.) 

According to section [3~T1 the ring C(s)[6q] is a Bezout domain. Consequently, there are elements 

a, b <G CIs, Sol such that 




(15) c = ap + bqe C[s, Sq] 

is a g.c.d. in C(s)[6q]. Hence, p/c and q/c belong to C(s)[6q]. In particular, there are n, € TZq, 
di € C[s] with gcd TC (n,, di) = 1 (i = 1, 2) such that p/c = ni/di and q/c = n-i/di. It follows pd\ = 
cni, qdi — cn-i. Consequently, both d\ and d 2 divide c in TZq. Since d\ and d 2 are polynomials, 
they possess a least common multiple h = di<i2/gcd(<ii, d 2 ), and it follows c = c/h £ TZq. Clearly, 
c divides both, p and q. Dividing (115|) by c yields the equation 

(16) a?iid 2 / gcd(di,d 2 )+bn 2 di/ gcd(di,d 2 ) = did 2 / gcd(di , d 2 ) . 

" v ' " v ' V v ' 

p q h 

By the coprimeness of n\ and d\, resp. n 2 and d 2l it follows gcd(p, h) = d 2 / gcd{d\,d 2 ), resp. 
gcd{q,h) — di/ gcd(di,d 2 ). Thus, gcd(q,h) and gcd(p,h) are coprime and, since by equation (fT6|) 
any common divisor of ji and q divides h, we can finally conclude the coprimeness of p and q. 
Thus, by Lemma [3~3l there are a*, b* £ TZq such that a*p + b*q = 1. The claim follows directly by 
multiplying this equation by c. □ 
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4. Controllability analysis 

4.1. Systems and Dynamics. 

Notation 4.1. Unless otherwise stated, the submodule spanned by a subset S of an R-module M 
is written [S] . 

Definition 4.1. An R -system A, or a system over R, is an R-module. 

Definition 4.2. A presentation matrix of a finitely presented R-system £ is a matrix P such that 
£ = [w]/[Pw] where [v] is free with basis v. 

Definition 4.3. An i?-dynamics, or a dynamics over R, is an R-system A equipped with an input, 
i.e., a subset u of A which may be empty, such that the quotient R-module A/[it] is torsion. The 
input u is independent if the R-module [u] is free, with basis u. 

Definition 4.4. An output y is a subset, which may be empty, of A. An input-output i?-system, 
or an input-output system over R, is an R-dynamics equipped with an output. 

Definition 4.5. Let A be an R-algebra and A be an R-system. The A-module A <Sir A is an 
A- system, which extends A. 

4.2. System controllabilities. In this section we emphasize several controllability notions which 
are defined directly on the basis of the above system definition without referring to a solution 
space. For the latter we refer to the next subsection. Let us start with some purely algebraic 
definitions: 

Definition 4.6 (see [13]). Let A be an R algebra. An R-system A is said to be A-torsion free 
controllable (resp. A-projective controllable, A-free controllable,) if the A-module A®nK is torsion 
free (resp. projective, free). An i?-torsion free (resp. i?-projective, R-hee) controllable R-system 
is simply called torsion free (resp. projective, freej controllable. 

Elementary homological algebra (see, e.g., [35]) yields 

Proposition 4.1. A- free (resp. A-projective) controllability implies A-projective (resp. A-torsion 
free) controllability. 

Proposition 4.2. R-free controllability implies A-free controllability for any R-algebra A. More 
generally, given any R-system £ that is a direct sum of a torsion module t£ and a free module A, 
the extended system A (&r £ is a direct sum of the torsion module A <&r t£ and the free module 

A(g> R A. 

The importance of the notions of torsion free and free controllability is intuitively clear: While 
the first one refers to the absence of a nontrivial subsystem which is governed by an autonomous 
system of equations, the latter refers to the possibility to freely express all system variables in 
terms of a basis of the system module. For this reason, and, secondarily, in reminiscence to the 
theory of nonlinear finite dimensional systems, we have the following: 

Definition 4.7. Take an A-free controllable R-system A with a finite output y. This output is 
said to be A-flat, or A-basic, if y is a basis of A®r A. Lf A = R then y is simply called flat, or 
basic. 

In finite dimensional linear systems theory, the so called Hautus criterion is a quite popular tool 
for checking controllability. This criterion has been generalized to delay systems (see, e.g., [30]) 
and to the more general convolutional systems defined over £' [42] and ftA$ [46]. All those rings 
may be embedded into the ring of entire functions via the Laplace transform. This motivates the 
following quite general definition: 

Definition and Proposition 4.1. Let R be any ring that is isomorphic to a subring of the ring 
O of entire functions with pointwise defined multiplication. Denote the embedding R — > O by Jzf '. 
A finitely presented R-system with presentation matrix P is said to be spectrally controllable if one 
of the following equivalent conditions holds: 

((i)) The O-matrix P = £f(P) satisfies 3k £ N : Ver £ C : trP(<r) = k. 
((ii)) The module £q = O ®r £ is free. 
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Proof. The ring O is an elementary divisor domain [19]. As a consequence, over O, any matrix 
admits a Smith normal form by left and right multiplication with unimodular matrices. Since the 
units in O are just the functions which possess no zeros in C, the rank of the Smith normal form 
equals the rank of P(c) for any <j € C Thus, the rank of P(a) remains constant if and only if the 
nonzero entries of the Smith normal form possess no zeros in C which, in turn, is equivalent to the 
absence of nontrivial torsion elements in [d]/[Pv], i.e., in Eq. □ 

Proposition 4.3. Let R be any Bezout domain that is isomorphic to a subring of O with the 
embedding R — ► O denoted by S£ . Then the notions of spectral controllability and R-torsion free 
controllability are equivalent if and only if _Sf maps non-units in R to non-units in O. 

Proof. Since R is a Bezout domain, torsion freeness of E implies freeness. Tensoring with the free 
module O yields another free module Eq, and, thus, by Definition and Proposition 14.11 spectral 
controllability. Again, since R is a Bezout domain, any presentation matrix admits a Hermite 
form. Thus, the torsion submodule tE of S can be presented by a triangular square matrix tP of 
full rank. If E is not torsion-free, at least one diagonal entry of this matrix is not a unit in R. If 
this entry is mapped to a non-unit in Eq by «£? , it admits a complex zero cxo- Thus, _£f (tP) has a 
loss off rank at a — ero- Contrary, if there is a non unit r € R which corresponds to a unit f £ O, 
consider E = [T]/[rr]. Obviously, the image of r in Eq is zero. Thus, the trivial module Ee> is 
torsion free. □ 

Remark 4.1. Note that, under the additional assumption that E admits a presentation matrix of 
full row-rank, the assumption of R being a Bezout domain may be replaced by a less restrictive one. 
In this case, equivalence of (Q <g>fl R) (1 O -torsion free controllability, with Q the ring of rational 
functions in one complex variable, and spectral controllabillity may be established (see, e.g., [30,46] 
for different examples) . 

We are now able to state the main result of our paper: 

Theorem 4.1. The convolutional system E defined in Definition \2.1\ is free if and only if it is 
torsion free. More generally E = tE © E/tE where tE is torsion and E/tE is free. Moreover, E is 
spectrally controllable if and only if it is torsion free. 

Proof. Recall that, according to Definition 12. 1[ E = 7?.<g>7?. Q SQ and 72<j is a Bezout domain by 
Proposition 13.11 Since the first assertion holds for finitely presented modules over any Bezout 
domain, it holds for Eq. The second assertion follows from Proposition 14.31 (The fact that the 
Laplace transform maps any non-unit of IZq to a non-unit in O is obvious.) Clearly, both results 
hold as well for E, which is obtained by an extension of scalars. □ 

4.3. Trajectorian controllability. In this section we will give two different interpretations of our 
algebraic controllability results that directly refer to trajectories of the system, i.e., to (generalized) 
functions which may be assigned to the system variables. To this end we need to introduce the 
notions of a solution space and a trajectory. 

Definition 4.8. Let E be an R-system and J? a space of generalized functions. The space & is 
called a solution space of E if it can be equipped with the structure of an R-module. 

Definition 4.9 (see [15]). Let & be a solution space of an R-system E. An J?" -trajectory o/E is 
an element o/HoniR(S, &). 

The crux of the first controllability notion (Def. I4.10|) is the possibility to assign an arbitrary 
(generalized) function from & to any system variable. 

Definition 4.10 (see [15]). An R-system is called ^"-trajectory controllable if for any element 
a £ E and any b £ & there exists a trajectory f with f(a) = b. 

The following result is borrowed from [15] and applies to any torsion-free controllable P-system 
where R C M. 

Proposition 4.4. The system E^/tE^ is M -trajectory- controllable. 

Another controllability notion is the following due to [45]. As the above it relies on the notion 
of a trajectory. However, since it refers to the possibility of connecting trajectories, the notions of 
future and past come into play. Thus, an appropriate solution space should allow the definition of 
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such local properties. This is not possible for the field of Mikusihski operators in general but for its 
subring .Mr, and more generally for the space B of Bochmians [1,26]. The controllability criterion 
in the behavioural framework is the possibility of concatenating trajectories. In our algebraic 
setting we may formulate this criterion as follows. 

Definition 4.11 (cf. [34,45]). Let E be an R-system and JP a solution space o/E that possesses the 
structure of a sheaf on R. Then E is called ^"-behavioral-controllable if for any two trajectories 
/i, fi £ Hom(E, JF) there exists f £ Hom(E, <&) such that for any a £ E there are £", t 2 £ R with 
/(o)l(-oo,tf) = /i(a)|(_oo,tf) and /(a)l(tg.oo) = /2(a)|( t g,oo)- 

Theorem 4.2. The system Er^Er, where Er is defined in Definition \2.1\. is B '-behavioural con- 
trollable. 

Proof. Since Er/^Er is free, any homomorphism is uniquely determined by the functions assigned 
to the basis. Thus, for the basis b = b\,...,b n we may chose t\ > t\ and set 

[J2{b), t > t° 2 

Moreover, any a £ Er^Er is given by a — X}j=n a ^ where the Ui have compact support. Thus 
there exist T\, T 2 such that suppo^ C [Ti,T 2 ], i = 1, . . . ,n. The claim follows by an application 
of the theorem of supports tf = t\+T u t% = t\ + T 2 (see. [1]). □ 



Remark 4.2. When distinguishing the cases a > and a = in ilb\) one could alternatively 
prove £ -behavioural controllability (resp. V -behavioural controllability) in the case a > or £ 2 - 
behavioural controllability (resp. T>' 2 -behavioural controllability) in the case a — 0, where £ is the 
space of infinitely differ entiable functions, T>' the space of Schwartz-Distributions, £2 the space of 
Gevrey- Functions of order less than 2, and T>' 2 the space of Gevrey ultradistributions. 

5. An example: two boundary coupled equations 

In order to illustrate our results, in the following we discuss a simple example. Consider the 
system of two second order equations 

(17a) d x Wi(x) = aw, i — 1,2, 

defined on an open neighbourhood f^ of [0, £i] C R, where a = as 2 + (3s + c. Those equations are 
coupled via the boundary conditions (i = 1, 2) 

(17b) HiiWi(£i) + ^w'Si) = 

(17c) w t (0)=u. 

According to Section W2\ the general solution of the initial value problems associated with ()17a|) 
reads (i — 1,2) 

, , fwi(x)\ _ f C(x-£i) S(x-ii)\ (cii 

[ ' WiWJ \crS(x-£i) C(x-U))\cn, 

with c,i = Wi(£i), ca = d x Wi(£i). The boundary conditions at x = £i yield 
(19a) HnCii + Mj2Cj2 = 

(19b) C(£ t )c a - S{£i)ci2 = u. 

Here, the relations S(—&i) — —S(£i) and C(—£i) = C(£i), derived in Section |2~^1 have already been 
incorporated. 

Thus, according to Definition 12.11 the convolutional system E associated with the boundary 
value problem (JTTJ) is the 1Z module [en, C12, C21, C22, w] the generators of which are subject to the 
equations (TH)]) . 

In order to reduce the number of equations, we aim to introduce new variables u>\ and u>2 such 
that (|19a|) is satisfied automatically, i.e., 

c-n = -Mi2W», c a = ^nuii, i = 1,2. 

Indeed, since 

W « = ~2 — I 2" ( - Mi2Cl +Mil c 2) ) i=l,2, 

Mil + Mi2 
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the new variables belong to S. Using the new generators wj., lo 2 and u, equation (|19b[) may be 
rewritten to obtain 

u = -^cjj, p t = iM2C(£i) + iHiS(£i), i = 1,2. 

Thus, piwi -j»2^2 = 0, and S = [lJi,lj 2 }/[pi^i — P2&2]. 

In accordance with Section f2.ll assume that Hi = Hit, with «; e N and i = 1,2. Thus, by 
Theorem 14. 11 checking spectral, torsion free, and free controllability are equivalent. Since the aim 
of this section is not the presentation of a general controllability analysis for the boundary value 
problem l|17|l but rather to give an example for the application of the derived algebraic results, we 
shall restrict ourselves to particular values for n\ and n 2 . In order to avoid tedious computations, 
we chose simply n\ = 1, n 2 = 2. Apart from that, we discuss the generic case only, i.e., we do not 
care about singularities which may occur for particular values of the /i.y, i,j = 1,2. 

Applying the algorithms of Section |3~T1 we obtain p\r\ + p 2 r 2 = e with 



r\ = 2((// 2 iA i ii - M22Mi20-)C(f) + (M22M11 - M21M12) <rS(£)) 

2M21M11M12 - M22M11 



r2 = n\%o - [i\ x 



e = -M22M12 (ct-ct), ct- 

M22M12 

Following Section 13.21 it remains to modify n, f2 in such a way that e is replaced by a constant. 
This may be done by applying the induction step of Lemma l3.3l once. To this end, let f 1 , f 2 , pi , P2 
be the complex numbers obtained by setting a = a in the Laplace transforms of ri,r2,Pi,P2- 
Assume that neither pi nor p2 are zero. Then the variables 

p 2 ri - fip 2 L g (pi)r 2 - L s {r 2 )pi 

1i — ) 92 — ~ 

P2£ Pie 

belong to IZq and, therefore, to 1Z. Thus, we have the Bezout equation piqi + p 2 q 2 = 1. 
From the above results, one easily verifies that with 

y = q 2 u>i + qiuj 2 

one has oj\ = p 2 y and uj 2 = p\y. Hence, y is a basis of the system under consideration. 



6. Conclusion 

For a class of convolutional systems associated with boundary coupled second order partial dif- 
ferential equations we have derived algebraic controllability results which translate directly into 
trajectory related controllability conditions. These results rely on a division algorithm for a partic- 
ular ring of Mikusihski operators with compact support that is obtained from the operator solution 
of the Cauchy problem associated with the given system of partial differential equations. However, 
this means that our algebraic setting does not apply directly to the given boundary value problem 
but rather to a convolutional system arising from these solutions in connection with the bound- 
ary conditions. A promising approach allowing an algebraic treatment from the very beginning is 
currently under investigation. 

The current work was motivated by previous contributions [2, 17] in which similar results where 
presented for differential delay systems. Those approaches have been shown to be useful not 
only for controllability analysis but also for the design of closed loop control schemes using the 
factorization approach or the method of finite spectrum assignment [3,4,18]. This suggests the 
investigation of similar methods for the class of systems considered within this contribution. 
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Appendix A. Representation of the operators S(x) and C(x) 

In this section we give interpretations of the operators S(x) and C(x) introduced in Section [2l 
Actually, we restrict to S(x) from which C(x) can be easily deduced by differentiation w.r.t. x. 
If a > in equation (|fb[) we may rewrite a as 



i ! i \2 i\ i— & b 2 c 

a = r (ia + a) -0), r = ^a, «=-, f 3 = ^ — - -. 

According to [23] the operator S(x) corresponds to the compactly supported function 



S(X, t) = I (h(t + XT) - h{t - XT))—J (P^T 2 X 2 -t 2 ) V , 

where Jo denotes the Bessel function of order zero. Thus, for any approximate identity (</> n )> the 
operator S(x) possesses a regular representation given by (cf. also App. IB.2|) 



(20) S(X) = f n (x)/<p n , fn(x,t)= J -—J ((3^T 2 X 2 - e)Vn{t ~ £)#■ 



XT 



Contrary, if a = in JTbJ), S(x) can be written as power series in the differentiation operator: 

oo 

S(x)=J2- 



(s - a) k x 2k+1 



t (2k + 1)1 • 

The convergence of this series is verified directly from the regular representation 

S ( x ) = 'V' 7K\ , i\i'" ' V'fc+l.n = 1pk,n ~ aipk,n, 1p0,n = ¥>n 

f~f (2fc+l)!v3„ 

where (ip n ) is any approximate identity of Gevrey order less than 2. 

Appendix B. Mikusinski Operators and Boehmians 

B.I. Generalized quotients. The set of locally integrable functions with left-bounded support 
forms a commutative ring £ + with respect to the pointwise addition and the convolution product. 
A celebrated theorem of Titchmarsh ( [40, Theorem VII], [41, Theorem 151] jj states the following: 

Theorem B.l. Assume that the convolution product of two locally integrable functions f and 
g the support of which is contained in M + vanishes on [0, T] . Then there exist nonnegative real 
numbers T g , Tf with T g + Tf ^ T such that both, f and g, vanish identically on [0, Tf\ and [0, T g ] 
respectively. 

Corollary B.l. The ring C+ is free of divisors of zero. 

Definition B.l. The field M. of Mikusinski operators is the quotient field of complex-valued locally 
integrable functions on R with left bounded supporQ. 

Definition B.2. Consider a commutative ring R together with an R-module M . Let IS. be a family 
of sequences of R such that: 

(1) // (tp n ), (ip n ) G A then {cp n tp n ) G A. 

(2) For f,g G M and (<p n ) G A the equality of sequences (f<p n ) = (gfn) implies f = g. 
Then the elements of A will be called A-sequences in M [25,28]. 

Definition B.3. Consider an R-module M with R a commutative ring and A a family of A- 
sequences in M . Let &/(M, A) the set of all pairs of sequences (f n ) G M N and (<p n ) G A satisfying 
(fifj — ipjfi for all i,j G N. With the (equivalence) -relation ~ defined by 

(((fn), (<Pn)) ~ ((ffn), (lpn))j ^ UpiQj = Ipj fi for all i, j G NJ 

Then the space B(M, A) of Boehmians on M is defined as £?(M, A)/ ~. For notational simplicity 
a Boehmian is simply denoted as /„/(/;„. The addition on B(M,A) may be defined according to 

fn/^n + 9n/lpn = [ipnfn + (fin9n) I ' [fn^n) [25, 28]. 



For several alternative proofs see also [22-24, 48] . 

Contrary to this definition Jvl is sometimes defined as the quotient field of the convolution ring of continuous 
functions with support in R+. However, in both cases the obtained fields of fractions are isomorphic. 
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In the following £(R) denotes the space of locally integrable functions on R and £o0&) (resp. 
£+(R)) the subset containing the functions with bounded (resp. left-bounded) support. With the 
pointwise addition and the convolution product £q (resp. £ + ) form a commutative ring and £ 
(resp. £+) is a £o- m °dule (resp. £+-module). 

Theorem and Definition B.l. [25, 26] Let A be the family of all sequences {f n ) in £o satisfying 

(1) 3C e R : Vn G N : f R \<p n (t)\dt ^ C 

(2) VneE:f R ip n (t)dt = l 

(3) Ve G R 3 n G N : Vn ^ n : supp y>n C [-e, e] . 
TVien we may define the following spaces of Boehmians: 

• TTie ring Mo °/ Mikusiriski operators with compact support is defined as B([Co]c , A) [lj. 

• The ring Mb. D A^o of regular Mikusiriski operators is defined to be B([£+]c + , A). 

• The elements of the space B([£]c ,A), which for short is denoted by B, are simply called 
Boehmians. 

Obviously, the spaces Mb, and B possess the structure of Mq modules. 

Remark B.l. The sets Mo and Mr regarded as commutative rings are clearly isomorphic to 
subrings of of M. Since the rings £+ and £q are free of divisors of zero, for these spaces the 
equivalence relation in \B.3\ could be replaced by: 

((/n), fan)) ~ ((5n)) (^n))J & Upi9j = ipjfi for some i,j G NJ . 

B.2. Divisibility in the ring of Mikusiriski operators with compact support. In this 
section we shall state some divisibility properties of the ring Mq. For proofs we refer to the cited 
literature. 

Proposition B.l. The mapping T a defined by the pointwise multiplication of a function f G 
£o with an exponential function t i— ► c at defines an isomorphism on £q which extends to an 
isomorphism on Mo [23,4-6]. 

Proposition B.2. The mapping L : £q — + C assigning to every element of £o the value of its 
integral can be shown to be a homomorphism. Its unique extension to Mo is denoted by the same 
symbol. 

Proposition B.3. An operator a G Mo divisible by (s + a), a G C if and only if L o T a (a) = 0. 

Remark B.2. Note that the function a : C — * C given by a(o~) = L o T^ cr (a) is the Laplace 
transform Jzf(a) of a. It is an entire function which satisfies a growth condition on the imaginary 
axis derived in [5]. Within this context, Proposition \B. 3\ means that s — a divides a in Mo if and 
only if a — a divides the Laplace transform of a within the space of entire functions. 
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